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ABSTRACT 

This paper describes an on-going simple laboratory experiment, referred to 
as the Beam Control Experiment (BCE) , which has the essential features of a 
large flexible structure. The experiment is used to develop and evaluate iden- 
tification and control algorithms which look promising in the active control of 
high performance large space structures. Some' results on the maximum likeli- 
hood identification of the parameters of the beam-actuator-sensor assembly from 
experimental data is presented in the paper. 

I. INTRODUCTION 

One of the major problems in the design of control systems which operate 
in the presence of a flexible structure is obtaining accurate information about 
the plant dynamics. In particular, knowledge of the frequencies, damping ratios 
and mode shapes c the flexible modes is critical to the successful design of 
a high performance system. System identification is an iterative process, the 
success ot which depends upon the choice of algorithm and system model, the 
choice of inputs to excite the system, and the quality of output data. A care- 
ful integration of these items is especially critical in the case of large 
flexible structures. 

In this paper we describe the development and performance testing of a 
simple laboratory model of a jitter control sys un designed to provide a stable 
image with optical components mounted on a flexible structure. The study will 
be carried out in three stages: (a) identification with simulated data, (b) 
identification with real data, and (c) comparison of closed loop performance 
with simulated results. Results from the first two stages are reported in 
this paper. 

This paper is organized as follows: Section 2 describes the experimental 

set-up and a mathematical model for the BCE is developed in Section 3. A brief 
description of the maximum likelihood estimation (MLE) algorithm is presented 
in Section 4. Results on the Identification of the parameters of the BCE using 
both simulated and experimental data are discussed in Section 5. A summary 
and future work is described in Section 6. 

II. BEAM CONTROL EXPERIMENT 

The idea behind this experiment (Fig. 1) is to demonstrate the inter- 
action between the control of an optical system, symbolized by a laser beam, 
and control of a flexible structure, represented by a flexible aluminum beam 
to which passive and active mirrors are attached. These mirrors bounce the 
laser beam toward a desired target. The inteiesting control problem stems 
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from the fact that the active mirror is in fact part of a proof-mass actuator. 
Thus, any attempt to control the laser beam will tend to disturb the aluminum 
beam, thereby also disturbing the laser beam. This intricate coupling is 
quite a challenge for a classical design but more amenable to modern techniques. 
The other aspect of the experiment is the use of a commercially available 
cated microprocessor (1SI MCP-100) capable of handling at maximum a 32- e 
Kalman filter at a 2000-Hz sampling rate. Such implementation, aside f 1 its 
laboratory usefulness, bring control technology one step further toward real 
space implementation, and the experience gained is valuable. 

The schematic of the experiment is shown in Figs. 2 and 3. The laser beau 
bounces first on a mirror situated near the middle point of the vertical alumi- 
num beam, then on a mirror attached to the tip of a pivoted proof-mass actua- 
tor (PPM). The laser beam finally reaches a photodetector, which measures the 
laser beam position. A beam splitter provides a visual display of the jitter 
by projecting the spot on a remote screen. 

Two sensors are used in the experiment: the photodetector measuring the 

liae-of-sight (LOS) error, and the TPM rate sensor measuring the relative velo- 
city of the proof-mass. 

A preliminary experiment had been performed earlier (Ref. 1) using com- 
mercially available software for identification and control synthesis. Only 
one sensor was used (position) and the control system was able to signifi- 
cantly damp out the beam vibrations; thus stabilizing the line-of-sight. 

However, in order to eliminate the static error and achieve a higher perform- 
ance, a better model is needed and thus more sophisticated identification 
techniques are sought to that end. 

For purposes of identification, a known control force is applied to the 
beam using the proof-mass actuator. The control force time-history and the 
beam position and relative rate outputs are recorded on a Nicolet 4094 digital 
oscilloscope. Special software transfers the input and output data from the 
Nicolet 4094 to a Harris 800 computer where the identification algorithms 
were exercized. Thus an efficient link between hardware ter'*s and sophisticated 
computer analyses (Tig. 4) was established. 

III. MATHEMATICAL MODEL 

In this Section, a state space model of the system is developed. The 
mathematical form of this model will be used both for simulation and identifi- 
cation of the parameters of the BCE. 

The angular displacement, 0 a , of the proof-mass actuator is limited to 
to a few degrees. For small angles the force and torque applied by the actua- 
tor are given by the equations (Reference 2) 


"a 

(i) 

mb# 

(2) 


a a 
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where m is the mass of the proof-mass actuator, b is the distance from the cen- 
ter of mass to the proof-mass pivot and I is the centroid inertia of the proof- 
mass actuator. The dynamics of tie aluminum beam will be defined in terms of 
the principal modes and mode shapes. Let be the modal amplitude of the i c ^ 

mode and define: 

= angular velocity of the i t ^ 1 mode 
th 

= damping of the i mode 

- translational mode shape of the i ^ mode at the 
beam tip (where the PPM is mounted) 

and <j>f = rotational mode shape of the i th mode at this tip. Let 4> im and $± m 
be the corresponding values for che mode shapes where the mirror is mounted 
on the flexible beam. The modax equations for the beam are: 

V i + 2 ^i w i + “i *i = ■ ^i f a “ *i T i = 1,2, .. M (3) 


where M is the number of modes represented in the model. Due to the actuator 
dynamics, the control force f applied to the actuator is related to f by the 
equation 


rf = rf - (K + mgb) (6 - 0 ) - D (0 - 0 ) (4) 

a c d l at 


where 


and 


r * lever arm of the actuator 
g * acceleration due to gravity 
K - spring constant of the actuator 
D * damping constant of the actuator 
0^ * rotation of the beam tip 


The rotation of the beam tip can be expressed in terms of the modal amplitudes 

by 


M R 

0 t = I \ q i ' 


(5) 
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Equations (1) - (5) can be reduced to the set of equations 

e* = - (K + mgb)/I e - D/I 6 + «C/I E q, + D/I E q, + 4 f (6) 

» « 3 XX 11 I C 


and 


■>1 ■ 


A i ®a - C i \ - 2 5 1 “i ’l - “l ’l + *1 E *i »i + c i £ *i ii 


+ B. f (7) 

ic v J 


where 



A, = 
1 

-K (I <f>* + mb((. i )/I 

(8) 


B. = 
1 

r/K 

(9) 

and 

C, = 
1 

A t D/K 

(10) 


Let yi be the displacement of the laser beam from the reference point \i,e., 
this is a measurement of the L.O.S. error). Let y 2 be the relative angular 
rate between the actuator and the beam tip. Then, 


y l ‘ 2 t ’i - (I 1 + l 2> 1 »lm + V. ] 

and 

y 2 = K - K m K - 1 q i 


(id 


( 12 ) 


where is the distance between the mirror on the beam and the mirror on the 
proof-mass actuator and ^2 *- s the distance between the photodetector and the 
mirror on the proof-mass actuator (See Fig. 3) . 

Equations 6, 7, 11, and 12 give a state space representation of the 
input/output behavior of the BCE with T0 a 0 a ^1 ^1 • • • ^mI as the 
state vector. , 

For a single mode model the equations are given by 


x * Fx + Gu (13) 

y * hx (14) 
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where 


and 


where 


Also, define 



0 

1 0 

0 

F = 

2 

~ui a 

-D/I a) 2 

z T 1 

D/I 


0 

0 0 

1 


" A 1 

- C 1 

c i*j - 2 'i“i 

G = 

1 1 

0 

r/I Ob/] 

f 

! 


h » 2 0 

2 ^lm ( *1 + V 

0 

L 0 i 

0 

"+1 i 


Ug = (K + mgb) /I . 


2/l(K-Hngb) 


(15) 


(16) 

(17) 

(18) 
09) 


£ 

The F, G and H matrices depend on the unknown parameters ( w a ,D,a>i,ei,4>i»4>]., 

♦ and the kn^wn parameters (l,m,b,r,£i and £2 ^ • The values of the 

known parameters are tabulated in Table 1* In this model the number of unknown 
parameters is equal to (6*H-2) where M is the numoer of modes. 


I 

4X1 0" 4 Kg-m 2 

m 

0.07895 Kg 

b 

0.06985 m 

r 

0.021 m 

*1 

0.165 m 

£„ 

0.217 m 

2 



Table 1. Known Parameters of the BCE 
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IV. MAXIMUM LIKELIHOOD ESTIMATION (MLE) 


There are several methods available for the estimation/ identification of 
parameters Figure 6 shows the main components of an identification method. 
The system to be identified and a mathematical model, M(p) , of the system 
is excited by a known input u. An error function, L(p,e), is formed from the 
outputs of the system and the model. Identification is the process of 
adjusting the model parameters p to minimize the error function. The choice 
of M(p) , L(p,e) and the adjusting mechanism foi p lead to different identifi- 
cation algorithms. In this paper, we shall restrict our attention to the 
maximum likelihood estimation of parameters. 

The MLE can be applied to a large class of problems and has good statis- 
tical convergence and accuracy properties. In addition, MLE is well suited 
for identifying physical parameters of the system. This is a drawback with 
most recursive algorithms. The main disadvantage of MLE is the amount of 
computation. However, the amount of computation can be reduced significantly 
by taking into account the special features of the dynamics of the large space 
structures . 

The flow of computation in the MLE is shown in Fig. 7. The mathematical 
model for the system is assumed to be 


x = F(p) x + G(p)u + (u 

(20) 

y = H(p) x + v 

(21) 


where x is the n-dimensional state vector, y the p-dimensional output vector 
and u the m-dimensional input vector, w and represent the process and 
measurement noise respectively. The matrices F, G and H depend on p, the vector 
of unknown system parameters. An input signal £u(t) , o < t < tjQ has been 
applied to the system and the output y of the system has been observed at dis- 
crete times t Q , tj, . . . t^. Further, it is assumed that x(o) and w are 
zero-mean gaussian random variables with 


cov 

(x(o)) 

= P(o) 

(22) 

cov 

(U>) 

- Q 

(23) 

and cov 

(V) 

- R 

(24) 

The identification problem 
experimental data u(t^) , 
y the output estimate, and 

consists 
y(ti) , 
e(t R ) be 

of estimating the parameters 
i *“ 1,2,...N. Let x be the 
the outout error 

p from the 
state estimate 


where = y ^k^ ~ Yv^)* ( 2 ^) 
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The negative log likelihood function, v(p,e), can be written as 

N t -1 

v (p,e) = - log L(p,e) = T. e(i) E(i) e(i) + log j B<i) j (26) 

i=l 

The maximum likelihood estimate of the parameters p is obtained by maxi- 
mizing the likelihood function L(p,e) (or by minimizing V(p,e)). This non- 
linear minimization problem has to be solved by numerical methods and makes the 
MLE computationally intensive. The computational aspects of MLE are discussed 
in Reference 4. 

V. IDENTIFICATION RESULTS 

Numerical results on the identit ication of the parameters will be pre- 
sented in two steps. First results from the identification of simulated data 
will be shown. This will be followed by results from experimental data. 

a. Simulation Results : Identif icaticn with simulated data was done 

to get a better understanding of the dynamics of the BCE, to pro- 
vide guidelines to set up the experiment and to test the MLE 
software. The system was simulated using 4 modes and was excited 
by a "bang-bang" type input with an amplitude of ±0.1 Newton. 

Figure 8 shows the laser beam position and relative velocity 
output from the simulation. This input/output simulated data 
was used to identify a single mode beam model of the system 

(see equations 13-17). The negative log likelihood function 
V(p,e) was probed at a few points to see its variation with param- 
eter p. Figure 9 shows the variation of the likelihood surface 
with w a and . The surface is well-behaved in these two variables. 
The damping terms £ a and were set to thg simulation values and 
only parameters (w a . and ^i m ) were allowed to 

vary. Table 2 shows the results based on simulated data. There 
is good agreement between simulated and estimated values. Now we 
are ready to try the identif ication with experimental data. 

b. Experimental Result s: The aluminum beam was excited by applying 

to the proof-mass actuator a sinusoidal force with a linearly 
varying frequency (so-called "chirp* 1 excitation) . Figure 10 shows 
the control force f c . The position and rate measurements are shown 
in Figure 11. As before, the single mode model will be used as a 
starting point for the identification of parameters (w a , D, 

$ 1 , <p . Initially, the MLE had convergence problems. 

These were related to one or more of the following causes: 1) large 

differences between the values of some of the actual BCE parameters 
and those of the original simulation, 2) bias in the input force 
and position measurements, and 3) error in rate measurement cali- 
bration. 
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The bias was accounted for simply by subtracting a constant from 
the force input and position output. The error in rate calibra- 
tion was taken into account by defining a scale factor a. This 
results in a new H matrix where 


2l 2 

0 

2( *i» - - < £ l + *2> *L> 

- 

0 

0 

a 

0 

R 


The scale facto^ was estimated along with the other 8 parame- 
ters. The estimated values are shown in Table 3,. The model 
parameters shown in Table 3 for "simulation values" were obtained 
from an extremely simplified model of the aluminum beam (canti- 
levered with a point mass at tip). Thus, it is not surprising 
that some of the values obtained from the identification are very 
different. In particular, values of rotational mode shapes are 
quite sensitive to local inertias and masses. 


H 

Simulated 

Data 

Estimate 

(Simulation 

Data) 

Estimate 

(Experimental 

Data) 

U) 

a 

46.36 

45.5 

37.18 

w i 

54.63 

54.58 

40.32 

C a 

0.005 

.005* 

85E-4 

5 1 

0.005 

.005* 

0.01 

♦l 

3.328 

3.49 

7.48 


-13.51 

-13.61 

-110.0 

♦i. 

1.068 

0.98 

21.0 


-10.6 

-10.3 

2.49 

a 

1.0 

1.0* 

0.80 





^These parameters were set to their simulation values. 


Table 3 Results with Experimental Data 
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figure 11s! 'vs a comparison between the measured c.nd estimated values of the 
outputs. The estimated values were generated for the pramcler set which 
resulted from identification using experimental data (Table 3). Figure 11a 
•hows the measured and estimated values of the position. Figure lib is a 
\ _ow-up of the same curve to show the difference between measured and estimated 
ilues. Figures 11c and lid show curves similar to 11a and lib for the rate 
easurement. There is very good agreement between the model output and the 
experimental data. 

This model will be validated by taking the direct approach. In the third 
stage of the experiment, a control system will be designed using the identi- 
fied model. The predicted behavior of the control system will be compared 
with its experimental behavior. 

VI. CONCLUSIONS AND SUMMARY 

In this paper, we have described a laboratory experiment which has the 
salient features of controlling an optical system located on a flexible struc- 
ture. The experiment will be used as a test bed for designing control and 
identification algorithms for large space structures. The parameters of a 
.Model suitable for designing a control system were identified using maximum 
likelihood estimation. The real test of a model is of course how well ic 
satisfies the goal of modelling. Currently, we are designing a control 
system based on this model and the results of this final stage will be 
reported in another paper. 
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FIG* 1 BEAM CONTROL EXPERIMENT 
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FIG. 5 BEAM EXPERIMENT GEOMETRY 
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FIG. 6 


COMPONENTS OF AN IDENTIFICATION SCHEME 



FIG. 7 COMPUTATIONS IN THE MLE 
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FIG, 10 BEAM EXCITATION INPUT 
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